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ABSTRACT 

The baryonic acoustic feature in galaxy clustering is a promising tool for con- 
straining the nature of the cosmic acceleration, through measurements of expansion 
rates H and angular diameter distances D\. Angle- averaged measurements of clus- 
tering yield constraints on the quantity D\/H. However, to break the degeneracy 
between these two parameters one must measure the anisotropic correlation function 
as a function of both line-of-sight (radial) and transverse separations. Here we investi- 
gate how to most effectively to do so, using analytic techniques and mock catalogues. 
In particular, we examine multipolc expansions of the correlation function as well as 
"clustering wedges" £(A/i, s), where fx = sm/s and s\\ is the radial component of sep- 
aration s. Both techniques allow strong constraints on H and D^, as expected. The 
radial wedges strongly depend on H and the transverse wedges are sensitive to Da- 
Analyses around the region of the acoustic peak constrain H ~ 20% better when us- 
ing the wedge statistics than when using the monopole-quadrupolc combination. How- 
ever, we show that the hexadecapole allows substantially stronger constraints than the 
monopole and quadrupole alone, as well as analyzing the full shape of £. Our findings 
here demonstrate that wedge statistics provide a practical alternative technique to 
multipoles, that should be useful to test systematics and will provide comparable or 
better constraints. Finally, we predict the constraints from galaxy clustering that will 
be possible with a completed version of the ongoing Baryonic Oscillation Spectroscopic 
Survey. 
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1 INTRODUCTION 

The clustering of matter is a powerful tool to probe the evo- 
lution of the Universe. In recent years, the baryonic acous- 
tic feature has been detected in the clustering of galaxies. 
Originating from pre-recombination plasma waves at red- 
shifts z > 1100, this feature is strongly detected in the tem- 
perature fluctuations of the Cosmic Microwave Background 
radiation (CMB). Its detection at low redshifts serves as 
an important confirmation of the cold dark matter (CDM) 
paradigm, and serves as a link between the late and the early 
Universe (|Peebles fc Yulll97Cl . iMeiksin et alj|l999h . 

The baryonic acoustic feature is also a practical cos- 
mic standard ruler llEisenstein fc Hulll998l . lEisenstein" et al.l 
ll998l . lE~isenstein et aljl999h . because the plasma waves left a 
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distinct imprint at the surface of last scattering (z* ~ 1100). 
In the CMB temperature fluctuations AT/T this sound 
horizon scale of ~ 150 physical kpc corresponds to ~ 1° 
in th e sky and is known t o an accuracy of < 1.5% (la 
level; iKomatsu etal] l2009h . Hence it can be used as a 
calibrated scale of the angular diameter distance Da(z*)- 
By measuring a related signature imprinted in the dis- 
tribution of matter at late times, one can perform geo- 
metrical measurements to determine angular diameter dis- 
tances Da at different z, as well as measure expansion 
rates H. The importance of this statistical tool is ampli- 
fied considering the recent discovery of t he acceleration of 
the expansion of the Universe (z < 1; iRiess et al ] 1 19981. 
iPerlmutter et all 1 19991 ). Although both Da and H are im- 
portant cosmic measurements, measuring H directly at var- 
ious z puts strong constraints on understanding the nature 
of the apparent acceleration of the observed Universe, e.g 
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through the so-called "dark energy" equation of s tate w 
(iBlake fc Glazebrookl |2003|, ISeo fc Eisensteinl [20031 . iLinderi 
12001 iGlazebrook fc Blakell2005l ). 

In the matter distribution, this signature appears as os- 
cillations at fc > 0.1/iMpc -1 in the power-spectrum, -P(fc), 
corresponding to a bump of excess overdensity in the two- 
point correlation function, £(s), at the characteristic comov- 
ing scale of ~ 100 /i -1 Mpc. This scale, although not equal 
to the sound-horizon, is closely related to it , with well- 
understood differences dMeiksin et al.l 1 19991. TSmith et al.l 
120081 . ICrocce fc Scoccimarrdl2008l . ISanchez et alj|2008h . ~ 



The baryonic feature is imprinted in s) and can 
be used to measure H and Da- For purely geometric rea- 
sons, radial clustering (i.e, clustering in the line-of-sight of 
the observer; fi = 1) contains H(z) information, where the 
transverse direction (n = 0) yields Da{z), where [i is the co- 
sine of the angle between the total separation of the galaxies 
and the radial direction. 

For S/N reasons, most studies have focused on the 
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wavelets l|Arnalte-Mur et al.l |2011| ) i n clustering of 



ies and galaxy-clusters of the SDSS (|York et al 
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Two Degree Field G alaxy Redshift Survey ( Colless et al.l 
120031 ). the WiggleZ (iDrinkwater et alj |2010| ) survey and 
the 6dF Galaxy Survey (| Jones et al.ll2009l ~ galaxy samples. 
The feature in the projected tw o-point function of photo- z 



samples has been analysed by (IPadmanabhan et al. 120071 . 
IBlake et al1l2007l . lEstrada et al.ll2009l . ICrocce et al.ll201ll l 

The baryonic acoustic feature in the monopole has been 
shown to constrain the cosmological information in combi- 
nation D\l H. The degeneracy of this combination limits 
its constraining power on expansion models. In this study 
we discuss techniques to use anisotropic clustering to break 
this degeneracy. 

Radial clustering measurements have been attempted 
on the Sloan Digital Sky Survey (SPSS) L uminous Red 
Galaxy sample (LRGs; iGaztanaga et al.ll2009l). as well as o n 
the much smaller volume MAIN sample ( Tian et al.ll20ldl ). 
Interestingly, both studies show strong clustering measure- 
ments, relative to the monopole and ACDM predictions, 
near where the baryonic a coustic feature is expected. How- 
ever, iKazin et al.l (|2010bh suggested a different interpreta- 
tion of these measurements indicating that the measurement 
could be the result of sample variance due to the limited vol- 
ume. 

A few studies have investigated using the informa- 
tion from the full P{n, fc) plane to co nstrain dark energy 
by usi ng geometric redshift distortions. lAlcock fc Paczvnskil 
i| 19791) describe how an intrinsically spherical system ap- 
pears anisotropic due to geometric distortions. They point 
out that by reconstructing the original spherical shape, the 
true cosmology can be obtained. This effect is manifested in 
clustering measurements, which are assumed to be isotropic 
(the cosmic principle). When converting the observed red- 
shifts to comoving distances, the observer is required to 
assume a fiducial cosmology. Choosing an incorrect cos- 
mology causes geometric redshift-distortions (in addition to 



the dynamical distortions, due to peculiar velocities of the 
galaxies). For this reason the observed clustering signal pro- 
vides an opportunity to apply the Alcock-Paczynski test and 
constrain the true underlyin g cosmology, assuming dynam- 
ical effects are understood (|Kaiserl Il987h . In practice, the 
baryonic acoustic feature plays an important role as it breaks 
degeneracies between amplitude uncertainties (erg, bias of 
tracer to underlying matter over-densities) and geometric 
shifts. 

iHu fc Haimar] (|2003h suggested disentangling H from 
Da by analyzing baryonic acoustic rings in the two- 
dimensional power spectrum. Focusing on phase shifts, they 
found that this technique, based only on geometric effects, 
can constrain the expansion rate of the Universe when ap- 
plied to galaxy and galaxy-clusters samples a t intermediate 
redshi fts z < 0.5, combined with CMB priors. IWagner et al.l 
(|2008f ) used mock catalogues at z = 1, 3 to demonstrate the 
usefulness of the technique, and show that light- cone effects 
do no t have a significant impact on the results. IShoii et al.l 
(|2009l ) argue that H and Da information is encoded in the 
full 2D shape, and present a generic algorithm that attempts 
to take into account dynamic distortions on all scales, as- 
suming non-linear effects are understood. For a first attempt 
to use the £(s _l,sh) plane of the SD SS LRGs to constrain 
cosmology, see lOkumura et al.l (|200Sl ). 

These studies, though, do not take into account the 
complexity of constructing a realistic, reliable (and invert- 
ible) covariance matrix for the full 2D measurement of the 
power spectrum or correlation function. Also, one short-term 
concern is the fact that near future surveys will have low S /N 
in the 2D plane. 

IPadmanabhan fc White! (|2008l ) investigated a more 
practical approach in which the 2D results are projected into 
one dimensional statistics. They proposed to break the H — 
Da degeneracy by combining AP analysis of the monopole 
of the correlation function, £0, with the quadrupole, £2- They 
argue that these measurements can constrain the warping 
in the 2D correlation function which is sensitive to Da ■ H, 
breaking the D\/H degeneracy obtained when probing di- 
lations in the monopole. 

Here we follow up on their analysis by: 

(i) Investigating the effects of higher order multipoles; 

(ii) Introducing a new alternative projection statistic in 
the form of clustering wedges; 

(iii) Performing the analysis, for the first time, in config- 
uration space (£(/i, s)). 

We propose to use wide clustering wedges £(A/x, s) to 
yield constra ints on H an d Da- A similar concept has been 
suggested bv lKazin et all (|2010bl ) (e.g., see their Figure 7). 
Tests performed here convincingly show that even a wide 
"transverse" wedge of < [i < 0.5 strongly depends on Da 
and a wide "radial" wedge of 0.5 < /1 < 1 constrains H. 
We show that these Afi =1/2 wedges do have intermixing 
terms, but these can be corrected for and are reduced with 
decreasing A/j,. 

The outline of the paper is as follows: in |5]we describe 
the statistical methods and the mock galaxy catalogues used 
in our analysis. In |3]we describe various theoretical aspects 
of z— distortions ( H3.ip . the implicatoins of the AP effect 
on basic cosmological parameters via determinations of H 
and Da ( £I3.2[) . the different clustering statistics used in our 
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analysis f £|3.3[) . and the extraction of H and Da informa- 
tion through the AP test ( H3.4p . In Sjjwe run two proof of 
concept tests. Using mock catalogues, we show that the ge- 
ometrically distorted signal can be retrieved from the true 
signal through the AP test f H4.ll) . and by doing so retrieve 
unbiased constraints on H and Da. We start with the ideal 
case in which all effects are known except for the AP effect 
in H4.2I and then, in ij4.3l we gradually add amplitude ef- 
fects. In H4.4I we investigate the uncertainties of H and Da 
as a function of the separation range used and compare the 
results obtained by the various ID projection combinations. 
In iJS]we present predictions for BOSS. Finally, fJS] presents 
a discussion, and ij7]our main conclusions 

Unless otherwise stated, we use the standard flat ACDM 
cosmology (Ok = 0, w = —1) with [Amo, f^bo 5 h, n 3 , as] = 
[0.25, 0.04, 0.7, 1, 0.8]. We test geometric distortions by vary- 
ing only the equation of state of dark energy w, when con- 
verting z to comoving distances. Our choices of distortions 
are explained in H3. 21 compared to expected degeneracies for 
various choices of curvature SIk, and matter density JImo- 
Unless stated otherwise, all distances hereon are comoving. 

To avoid semantic confusion, we briefly explain here 
the terminology of the different spaces we explore. First, all 
analyses are based on two-point correlation functions, which 
we refer to as configuration-space, as opposed to k— space. 
Second, because geometric redshift distortions and dynamic 
redshift distortions are different in nature, we minimize the 
use of the generic term for both, "redshift distortions" , and 
call them "geometric distortions" and "dynamical distor- 
tions". Because we analyse geometric distortions with and 
without dynamical effects, hereon we avoid using the com- 
mon expression redshift- space (z — space). Instead, when dy- 
namical effects are applied we refer to it as velocity-space, 
and when they are not we refer to it as real-space. 



2 METHODS OF ANALYSIS 



2.1 Statistical tools 



In our analyses we use the lLandv fc Szalavl (|l993l ') £ estima- 
tor. For details of usage please see Appendix [El 

When constraining parameters, we use the standard x 
technique, where 



x 2 ($) = y, ( M * (*) - D i) <V(- W *(*) - D > 



(i) 



where i,j are the bins tested. For reasons described below, 
the "data" D is given by the distorted measurement, , 
meaning £ measured from our mock catalogues when us- 
ing the incorrect cosmology to convert redshifts to comov- 
ing distances (the AP effect). The base-template used for 
modelling is £ T , meaning the actual true clustering signal. 
Hence, the models M are given by the shifted measurements 
£ 5 , that is the result of using the parameters $ to shift the 
template £ r . 

When limiting $ to [H, Da] we calculate \ 2 using brute 
force on a 2D grid. When investigating a larger parameter 
space, we apply a Monte-Carlo-Markov-Chain (MCMC). We 
verify that both methods yield similar results. 

The statist ics we use have cov ariant uncertainties, (e.g., 
see Figure 4 in iTaruva et al]|201ll for the correlation coeffi- 
cient between £o and £2-) For this reason we define D and 



M in array format. For example, when analyzing monopole, 
quadrupole [i = 0,2) the data is defined as D = [£o)£a] = 

We construct the covariance matrix dj from the 
A'mocks = 160 mock true signals. (For mock description see 
£|2.2H When using multipoles (or wedges) combination 
we define the covariance matrix as: 



a 



1-3 



mocks 



m,n— 1 



(2) 



2.2 Mock galaxy catalogs 

To simulate the observer's point of view, we analyse the 
mock galaxy catalogues from LasDamas and the Horizon 
Run. A similar version of these mock catalog ues has been 
used i n our previous analyses of the monopole ( Kazin et al.l 
bOlOal ) and radial clustering (|Kazin et al.ll2010bF 

The LasDamas simulations use a cosmology of 
[fijv/o,^60,n s ,/i,o- 8 ] = [0.25,0.04,l,0.7,0.8] and the Horizon 
Run uses [0.26,0.044,0.96,0.72,0.8], where Q M is the present 
baryonic density and n 3 is the spectral index. Both these 
cosmologies are well motivated by constraints obtained by 
WMAP 5- year measurements of te mperature fluctuation s in 
the cosmic microwave background (|Komatsu et al j2009l ).To 
understand effects of velocity-space we analyse all volumes 
both in velocity- and real-space. 

The LasDamas collaboration provides realistic LRG 
mock catalogue^] by placing galaxies inside dark mat- 
ter halos using a Halo Occupation Distribution (HOD; 
iBerlind fc Weinberg! [20021 ). HOD parameters were chosen 
to reproduce the observed number density as well as the 
projected two-point correlation function w p (r p ) of galax- 
ies in the SDSS-LRG sample at separations 0.3 < r v < 
30fe _1 Mpc, below the scales considered here. For more de- 
tails see McBride et al. (in prep.). We use a suit of 160 LRG 
volume-limited mock catalogues constructed from light cone 
samples with a mean number density of n ~ 10~ 4 /i 3 Mpc -3 . 
Each mock catalogue covers the redshift range 0.16 < z < 
0.44 and reproduces the SDSS angular mask, giving a vol- 
ume of 1.2Gpc 3 /i -3 . The LasDamas real-space catalogues 
are similar to the velocity-space catalogues in all aspects, 
with the exception of the shift in z due to peculiar veloci- 
ties. 

The Horizon Rur0 provides an ensemble of 32 BOSS 
volume realizations of mock LRG samples with a higher 
number density than DR7, n ~ 3 x 10 -4 /i 3 Mpc~ 3 , as ex- 
pected in BOSS. LRG positions are determined by identi- 
fying physically self-bound dark matter sub-halos that are 
not tidally disrup ted by larger structures. For full details see 
I Kim et al. I (120091 1. We construct these mock catalogues by 
dividing each of the eight full sky samples of the Horizon 
Run into four quadrants. We map real-space into velocity- 
space, and limit the samples to the expected volume-limited 
region of the BOSS LRGs (0.16 < z < 0.6). 



1 http:/ /lss. phy.vanderbilt.edu/lasdamas/ 

2 http:/ /astro. kias. re. kr/Horizon-Run/ 



© 0000 RAS, MNRAS 000, 000-000 



4 Kazin E., Sanchez A. & Blanton M. 




Figure 1. Mean two-dimensional correlation functions £(/^, s) from the ensemble of mock catalogues in real- (left) and velocity-space 
(right). The contours are spaced logarithmically between 5 and 10 — 4 for positive values, and at —0.001, —0.005, —0.01 for negative values. 
The contours corresponding to the values 0, 1 are colored with thick yellow lines. The solid contour lines, following the color scheme, 
correspond to the result obtained when using the correct cosmology when converting z to comoving distances. The dashed lines show 
the geometrical distortions obtained by assuming = —1.1 instead of the true value uP~ = —1. It can be clearly seen that dynamical 
effects dominate over the geometric. 



3 THEORY 

3.1 Redshift distortions: geometric vs. dynamic 

Redshift distortions arise due to two effects when converting 
the redshift z h s of a galaxy into a comoving distance: 



X = c 



H(zy 



(3) 



The first effect involves the assumption that the ob- 
served redshift is produced entirely by the expansion of the 
Universe z cos . This assumption is, of course, incorrect in 
the presence of peculiar velocities, which introduce an ad- 
ditional Doppler component z pcc leading to radial shifts in 
the infered distances. Although these shifts are small com- 
pared to the true distance x{ z cos) (less than 1% at z ~ 0.3), 
they strongly affect clustering measurements which depend 
on separations between galaxies. We refer to these as dy- 
namical distortions. 

Another, more subtle, redshift distortion effect arises 
due to the conversion of redshift to distance using only ap- 
proximately known cosmological pa rameters. The c onver- 
sion relies on the Hubble parameter |Friedmanlll923 ) 



H(zf = Hi 



n M0 (i + zf + n K (i + zf 



±0, 



(4) 



where ili are the standard cosmological density terms 
at present day for matter (M0), curvature (K) and 
dark energy (DE). The Hubble constant Ho = H(0) 
jHubble fc HumasorJ[i~93ll ) factors out trivially and we thus 
express comoving distance in units of h~ Mpc, where h = 



_ffo/(100 kms _1 Mpc _1 ). The rest of the parameters have 
more important, and potentially measurable, effects. We re- 
fer to these AP effects as geometric distortions. 

One way of overcoming these effects is to recalculate 
clustering statistics for every set of parameters when de- 
termining cosmological constraints. However, that approach 
is currently not practical. Instead, we calculate £ using a 
fixed fiducial set of parameters, and vary the result using 
linear equations. As we show below, this method is accurate 
enough. 

Figure [T] illustrates dynamic and geometric distortions 
in the LasDamas mock catalogues using the anisotropic £ in 
the fj, — s plane. The information in this coordinate choice 
is similar to that in the commonly used sn — s± plane. 
We define s to be the spatial separation vector with ra- 
dial and transverse components su, s±. In real-space (left 
panel) the true signal corresponds to flat horizontal contour 
levels in £(/x,s), shown as colored contours (solid lines) A 
noticeable signature is the baryonic acoustic feature around 
s ~ 110/i _1 Mpc. 

The dashed lines show the result when we introduce 
geometric distortions by using w — —1.1 instead of the true 
value w = — 1 in converting redshifts to comoving distance. 
These distortions are more noticeable at large scales, though 
they are also present on small scales. 

The right-hand panel illustrates the equivalent measure- 
ments with the addition of dynamical distortions (velocity- 
space). It can be clearly seen that the dynamical distortions 
dominate over the geometric ones. 

Three noticeable features are worth mentioning here. 
First, the velocity dispersion effect is clearly seen in the 
clustering signal along the line of sight (fj, = 0). Although 
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commonly regarded as a small scale eff ect, it is still presen t 
on scales of 60 h~ 1 Mpc, as discussed bv lScoccimarrol (|2004h . 

Second,the negative sea along the radial direction is ap- 
parent at s ~ 60/i _1 Mpc in this cosmology. Notice that in 
real-space (left plot) £ turns negative only at ~ 135 h~ 1 Mpc. 

Third, the baryonic acoustic feature, which appears as 
a positive stripe in real-space clustering, appears here as 
ridges which decrease strongly in amplitude towards the line- 
of-sight. Here the radial baryonic acoustic peak is negative, 
but can be positive depending on the value of the squashing 
par ameter p = f l b. 

iKaiserl (| 19871 ) originally describes linear dynamical dis- 
tortions by coupling the logarithmic rate of change of the 
growth of structure / to /i. By doing so he related the under- 
lying real-space isotropic P(k) to the apparent anisotropic 
velocity-space one. The bias b is introduced when relating 
to matter tracers. In this study we focus on the more sub- 
tle ge ometric distortions, and refer the reader to iHamiltonl 
(199§) for a review of dynamical distortions. 

3.2 The cosmological power of the AP effect 

Throughout this study we explore techniques to break the 
geometric H — Da degeneracy with clustering. Here we ex- 
amine how these constraints are related to fundamental cos- 
mological par ameter s assu ming a ACDM model. 

Following iHoggl (|l999T ) we define the transverse comov- 
mg distance as: 



D 



sinh I 



X 

D H - 



for f^K > 0; 

for Q K = 0; (5) 

for Q K < 0, 



where Dh is the Hubble distance c/Ho, and is related to the 
angular diameter distance by Dm = Da(1 + z). As we ex- 
plain in £13.41 the AP effect can be quantified by the dilation 
and warping parameters a and e, which depend, in turn, on 
both H and Da (see Equations 1151 [161) . These quantities 
depend in a non-straightforward fashion on the density pa- 
rameters Qi, and the dark energy equation of state w. H 
is given by Equation Q, while Da depends on H and Q,k 
according to Equations ^ and © . 

Figure [2] shows how H and Da depend on cosmological 
parameters for a number of redshifts. In each panel (each 
redshift) we hold two of the three parameters JImo, ^k, and 
w (where Qa = 1 — £Imo — ^k) fixed to a "true" value and 
modify the third from its fiducial according to the fraction 
indicated in the legend, between 1 to 1.5. We clearly see that 
at low redshifts H and Da yield degenerate constraints on 
w, J7k and Qm that can be broken as z increases. We notice 
that the dependence on Qmo does not vary much as a func- 
tion of redshift, where both J1k and w align with the Da 
axis at high z, meaning H is not sensitive to these parame- 
ters. 

This plot demonstrates that the fii, w degeneracy can 
be broken when applying the AP effect at high redshift (z > 
2). 

In this study we examine AP effects when varying w at 
mock mean redshifts (z) = 0.33, and 0.44 as indicated by 
thick boxes in bottom panels of Figure [2] Figure [2] demon- 
strates that our results are similar to those we would have 
obtained by choosing to vary fix or JImo- 



3.3 One-dimensional projections of £(/t, s): 
introducing clustering wedges 

We define clustering wedges as 



(6) 



where fi is the cosine of the angle between the total separa- 
tion of the galaxies s and the line of sight. We assume here 
the plane-parallel, or small angle, approximation according 
to which two galaxies at the same distance from the observer 
yield /i = irrespective of their angular distance. We note 
that the baryonic acoustic feature scale at z = 0.3 corre- 
sponds to ~ 7° i n the sky, and is smaller at larger redshifts. 
ISamushia et al.l (|201ll ) discuss observer angle effects. Due to 
our methods of building templates, our models incorporate 
large angle effects, and we do not test for them. 

Using spherical harmonics, the anisotropic £(/x, s) may 
be written as: 

£(/i,«)= Y, PtMtt(*), (7) 

even I 

where Ve are Ledgendre polynomials (e.g, Vo = 1, Vi = 
|(^ 2 - 1), Ta = |(35^ 2 - 30/1 + 3)) and 

2/ + 1 r 1 

6 = ^^ / Vl(Mfa,8)dn. 



(8) 



Equations and can be used to find the relation 
between the clustering wedges and the multipoles. Discard- 
ing contributions from multipoles with I > 2 this relation is 
given by: 



f(A/x,s)=£o- 



1 ( Mmax Mm: 



16 



(0) 



A hexadecapole term would mean an additional term 
given by 



1 (7 ((& 



8 



Mmiri 



10 (/t 



3 - u 3 ■ 

max Pmin 



/'■■ 



/'■■ 



+ 3 & (i°) 



on the right hand side of Equation @, and higher multipoles 
can be calculated in a similar manner. 

For simplicity, in this study we focus on clustering 
wedges defined by a width of A/i = 1/2. Of course, this anal- 
ysis can be generalized to various wedge widths. We discuss 
the results obtained with various values of A/i in Appendix 

El 

Defining the radial wedge £|| as that given by 0.5 < /t < 
1 and the transverse (i as 0< fi < 0.5, Equation © yields: 



3 

Ja 



6 



(11) 



(12) 



The hexadecapole term would add a third column in 
the matrix on the right side of Equation (|11[) with absolute 
values of 15/128 ~ 0.12. 

If £(/t, s) consisted only of t = 0, 2 terms, the two 
A/j, — 0.5 wedges would form a complementary basis to 
that of the multipoles. In the more generic case, these wide 
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Figure 2. Each panel displays the degeneracies between $=[ui, t^MOi Q~k) in the H — Da plane as a function of z. The axes, are in 
units of a "true" cosmology <E> = [—1, 0.25,0]. The "false" values are given in fractions in increments of A(frac) = 0.005. The legend 
reads such that: w < means w = w T ■ frac, w > w T means w = I frac, f2jyio > ^mo means ^M0 = ^mo ' / mc > ^M0 < ^mo 
means Cmo = ^Mo// rac > > means frac — 1 and fix < means 1 — frac. In some high z panels we highlight 5, 10% deviations 
in <3>. The thick boxes at low z indicate the geometric effects we test in this study. The figure clearly shows that at (cosmologically) low 
redshifts there is a large degeneracy between the parameters. This is relaxed at higher redshifts where distortions in S7k and w effect 
mostly Da, and distortions in Qmo yield similar results to low z. 



clustering wedges comprise an alternative, but not totally 
complementary basis. It is easy to see that given any com- 
bination of even is, the monopole is always the average of 
the A/x = 0.5 wedges, but the quadrupole is combined with 
higher order multipole terms in a complicated fashion. This 
means that given non-zero £f>2 terms, these wide wedges 
do not contain exactly the same information as [{0,62], and 
hence form an alternative, non-complementary basis. To 
have a fully complementary basis to £ which contains N mul- 
tipoles would, of course, require the same number of wedges 
(or any other projection). 



In Appendix [D] we test the relationships between the 
clustering wedges and multipoles. We find that the two wide 
clustering wedges (Afi = 0.5) are defined fairly well by the 
monopole and quadrupole in velocity-space (and monopole 
only in real-space), and hence may be used as an alternative 
basis to these multipoles to project most of the information 
contained in £(/i, s). In the next section we utilize this fact to 
show the effectiveness of the wedges to understand geometric 
distortions, and use them to constrain H and Da- 
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3.4 Dilation and warping in clustering: a 
treatment of multipoles and wedges 

Here we show that radial clustering wedges are, as expected, 
mostly sensitive to H while the transverse ones are most 
sensitive to Da, even for two wide A/i = 0.5 clustering 
wedges. 

iPadmanabhan fc White! (|2008l ) parameterize geometric 
distortions in clustering. We make use of their Equations 
(2) — (4), and introduce them here in configuration space. 
We define s to be the true spatial separation vector with 
radial and transverse components s\\, s±. The geometrically 
distorted separatio ns are indicated by a T> s u persc ript. 

As shown by IPadmanabhan fc White] (120081 ). distor- 
tions to the components of the separation can be param- 
eterized by a factor a which causes isotropic dilation and a 
parameter e that causes anisotropic warping, such that: 



S||«(l + e)' 
s±a(l + e) 



(13) 
(14) 



The Jacobian of transformation between the true vol- 
ume element d 3 s and the distorted d 3 s E ' is o? . Given that 
the comoving separation dx = cdz/H(z), and that the phys- 
ical angular diameter distance is (1 + z) Da = X0 it is easy 
to show that the dilation parameter is given by 



H_ 



23 \ 1/3 



Da 
U A 



2/3 



(15) 



Applying Equation (fl5|) to Equations lfl3)) and (fl4)) 

yields 



l + e = 



H V D V^ 1/3 



HDa 



(16) 



The combination of Equations (|13[) and (j 141) yields: 
3 d = a (1 + 2eV 2 (m)) s, 



(H ) 2 = fi 2 + 6e(it 2 - //). 



(17) 
(18) 



Note the difference in signs betwe en configuration space \x 
and fc— space Hk (Equation (3) of IPadmanabhan fc White! 
120081 ). 

Substituting these last two equations into Equation ([7} 
yields: 



-, 6 \ , , , 4 d&(s) 



2r 



dln(s)' 



(19) 



(20) 



Here we neglect terms of order 0(e 2 ). See Appendix [Al for 

inclusio n of hexadecapole terms. 

As IPadmanabhan fc White! (|2008l ) mention, the second 
and third terms on the right hand side of Equation (|19[) 
effectively cance l each other out, leaving (. n '(s) ~ £o(&s). 
lEisenstein et alj (120051 ) and I Sanchez et al.l |2009l ) demon- 
strated that this relationship works very well on the SDSS 
DR3 and DR6 LRG samples respectively, showing that the 



monopole alone constrains the degenerate combination in a, 
meaning D\/H. 

IPadmanabhan fc White! |2008l ) showed that the com- 
bined information of £o(s) and £2(3) can be used to measure 
simultaneously a and e. Because these parameters depend on 
different combinations of H and Da, this can in turn break 
the degeneracy between these parameters obtained from an 
analysis based only on the monopole. Here we present a sim- 
ilar concept based on clustering wedges. 

By combining Equations JTT} with Equations (|19[1 and 
(|20|) it is possible to quantify the effect of geometrical dis- 
tortions on the clustering wedges: 



H 1 



£(») = a +Cx(c), 



(21) 
(22) 



where we have used the fact that for small e, a(l + 2e) « 
H v ' /H and a(l — e) ~ Da/ Da- These equations hold for 
clustering wedges in general, where for A/j, = 0.5 the correc- 
tion terms are given by 



= e - 



140dln(s) 140 ? 1 ' 
373 d£±(. 



5 d£ (8) 
4 dln(s) 
677 rfg||(s) _ 
840rfln(s) 840 dln(s) 

(^)) 



and 

Cx(e) 



35 



1 dCo(s) 53 dg 2 (s) 123 
4dln(s) 280 dln(s) 140* 2laSj 



(23) 



107dg||(s) 317dg ± (s) 
840 dln(s) 840 dln(s) 



+ 



..|fe,< 



(24) 



assuming 



flatness, sec i|3.2l for a more generic treatment 



We neglect Q(e 2 ) a nd higher contributions. Equation 25 of 
iTaruva et all (|201Ch gives a more generic treatment of the 
linear AP effect in the P(fj,, fc) plane, whereas the equations 
presented here are ID projections. 

In £14.1 1 we show the validity of the equations presented 
here. In i)4.2l we demonstrate that correcting for the AP 
effect yields H and Da to high accuracy, and compare the 
wedge technique to multipoles. 



4 PROJECTIONS IN PRACTICE: TESTING 
THE AP EFFECT 

Here we demonstrate the applicability of Equations (|21|l 
and (1221) using analytic formulae and mock galaxy cata- 
logues. We show that, as expected from these equations 
the wide A/i = 0.5 "radial" clustering wedge dominantly 
constrains H, and the wide "transverse" one is sensitive to 
Da- This means that the information from these clustering 
wedges breaks the degeneracy in the combination D\/H 
obtained from the monopole only, and that the underly- 
ing values of H and Da can be obtained at high accuracy. 
Here we also compare the results obtained in this way with 
those recovered from the alter native multipole technique of 
IPadmanabhan fc White] (|200Sl ). 
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As described in £|2. II we define the true £ r signal to be 
the mock mean results obtained using the true simulation 
cosmology when converting redshifts to comoving distances 
X- The distored signal £ 73 is similar to £ T except that we use 
a different cosmology to convert z into X- We perform this 
AP effect both in real- and velocity-space, and thus we apply 
geometrical distortions in both cases. Finally, we define the 
shifted signal £ s to be our attempts to reconstruct £ 23 from 
£ r . Technically this means that we transform £ T to £ s 



using Equations (|21[) and ([22} for the wedges and Equations 
(fT9|) . (f2(?J> for multipoles. 

The tests we perform are as following: 

(i) EH1 We show a near perfect dependence of the radial 
wedge on H and the transverse wedge on Da by shifting £ T 
results to match a £ x> signal. 

(ii) H4.2I We use H and Da as varying parameters when 
fitting a model constructed from a template (the £ T signal) 
to match "data points" (the £ 23 signal) , and show that the 
best fit constraints on these parameters agree with the true 
values. 

Our mechanism is s imilar in concept to that used by 
IPadmanabhan fc White! (|2008l ). with the difference that we 
simulate the observer's point of view by including large- 
angle effects and explicitly use a wrong cosmology when 
converting redshifts to c omov ing distances (Equation [3]). 
IPadmanabhan fc White! (|2008f ) warped distant boxes ac- 
cording to a given value of e. One main difference is that we 
focus on low values of the warping parameter because our 
derivations are valid for small e. Here we focus on results ob- 
tained by changing the dark energy equation of state from 
its true value w = — 1 value to —1.1 (yielding a = 0.9832, 
e = —0.0033 at the mean redshifts of the mocks, (z) — 0.33) 
as well as —0.9 (yielding a — 1.0175, e = +0.0035) to exam- 
ine two directions of shift. These choic es are semi-arbitrary , 
as w is known to an accuracy of 10% (iKomatsu et aLll2009l. 
ISanchez et aT1l2009l . iPercival et aHl20ld , iReid et alj|201ol ). 
The distorted cosmologies analysed here correspond to the 
squares shown in Figure in the z = 0.33 panel. For this 
low redshift these variations are highly degenerate with mis- 
estimating Qmo or Qk ■ Lastly, the radial direction used here 
(n = 1) is the bisecting vector of s originated from the ob- 
server. 



4.1 Analyzing geometric distortions: proof of 
concept using mock catalogs 

In this section we test the accuracy of Equations (|19[) - l|24p 
using the suit of 160 SDSS-II mock galaxy catalogues de- 
scribed in 3231 

We measure the correlation functions using the true cos- 
mology of the simulations (£ r measurements) and the in- 
correct value of w = — 1.1 (£' D measurements), and we use 
these equations to shift the £ r measurements to match the 
£° ones (£ s measurements). 

Our results are shown in Figure [3] The upper plots cor- 
respond to the results for the clustering multipoles and the 
bottom ones to the clustering wedges. The left plots show the 
measurements in real-space, and the right in velocity-space. 
Each plot consists of two panels. The top panels show £ • s 2 
The £ T results are shown by solid black lines, the £ v re- 
sults (called AP) are dashed red lines, and the £ s results 



are the dot-dashed lines. In the bottom panels, the £ D re- 
sults form the reference to which we compare differences of 
£ r (black) and £ 5 (blue) in units of the uncertainty . In 
the multipoles the monopole results are in diamonds, and 
the quadrupole results are squares. In the wedges the radial 
wedge (Afi > 0.5) results are in diamonds, and the trans- 
verse wedge (Afi < 0.5) results are squares. 

In the top left plot we see the AP effect on the multi- 
poles in real-space. We verify that the shift in the monopole 
from the £o~ signal (solid) to the f(P (dashed) is described 
very well at zeroth order in e, meaning by fio'(s) ~ f c | (as). 
Adding the warping e terms adds little. In velocity-space, 
however, we do notice improvements when adding the first 
order correction at s < 50/i -1 Mpc. 

In real-space we do not expect signal in higher order 
multipoles. In Figure |3j though, there is a slight detection 
of £2 and even (,e>2 measurements. These probably arise due 
to Poisson shot noise (either in the random points or data), 
or due to large angle effects. Nevertheless, we see that the 
AP effect is understood for £2 (Equation I20[) in both real- 
and velocity-space. As expected in £2, the a shift only is 
not sufficient, and adding the e terms (blue) explain the AP 
effect to high accuracy. We also test higher order corrections 
of £4, dln(£ s )/ds in the quadrupole and find them negligible. 

As for the £4 statistics, we find that in velocity-space 
our corrections work well at s > 30/i -1 Mpc. We argue that 
it does not work at smaller scales because we do not use the 
expected £e>4 terms, which are required due to leakage of 
multipoles in the AP effect. 

We find similar trends for the clustering wedges statis- 
tics (bottom plots). In real-space (bottom left) we expect 
both £ T wedges to coincide. We notice minute differences 
at 130 /i _1 Mpc, (amplified in the plot by s 2 which make 
them visible). In velocity-space, as expected from the large 
quadrupole, there is a clear separation between the wedges, 
where the radial wedge is strongly suppressed and the trans- 
verse wedge elevated. The blue symbols indicate that the AP 
is very well described by Equations (|21|l - (|24| l. The first order 
correction terms C| | : j_ (e) (lines; see legend) are small (< 1% 
on most scales) with respect to the £ r wedge signals. This 
indicates that even with our definition of wide "radial" and 
"transverse" wedges (Afj, = 0.5), the radial wedge is mostly 
sensitive to H, while the transverse wedge to Da - That said, 
in Appendix [Bl we show that although setting C\\ t ±(e) = 
yields fairly accurate H, Da results, including C| | ,_l (e) can 
improve the results. 

To summarize this test, we have proven here that AP 
effects on the clustering multipoles and wedges are very well 
understoo d. We show that the dilation and warping terms 
shown by IPadmanabhan fc White! (|2008f ) explain this effect 
very well in the monopole and quadrupole of the two-point 
correlation function. In addition we show an alternative ap- 
proach to measuring the AP effect by analyzing clustering 
wedges. We show here that even a wide "radial" wedge is, 
as expected, mostly sensitive to H and a wide "transverse" 
wedge to Da - In Appendix[C]we perform similar tests on an- 
alytic formulae and obtain similar conclusions. In the next 
section we investigate the power of this method to obtain 
constraints on H and Da from measurements of the multi- 
poles and wedges. 
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s [h _1 Mpc] s [h^Mpc] 

Figure 3. Mock mean clustering projections (based 160 LasDamas realizations) with and without the AP effect. Left: real-space. Right: 
velocity-space. Top: multipoles (applying Equations 1191 120)1 . Bottom: A/i = 1/2 wedges (Equations I2H24) |. In the top panel of each plot 
are the ID projections, and in the bottom panels the difference of each result with the DISTORTED signal (labeled as AP) normalized by 
the uncertainty of one realization. As indicated in legend, black solid lines/symbols are the true projection signals (T), the long-dashed 
red lines are the AP signals (£>). We apply the AP shift to the true signals to obtain the triple-dot-dashcd blue lines/symbols (5). A 
perfect shift would yield a null result for the blue symbols in the bottom panels. The dot-dashed purple and dashed orange lines, as 
indicated in the legend, are the first order e correction terms. The AP distortion applied here is using ui FID = —1.1 instead of the true 
value —1 when converting z to comoving distances. 



4.2 Reproducing the true H and Da 

Here we perform the AP test on the LRG SDSS-II mock 
catalogues described in £12.21 to measure H and Da- When 
quoting uncertainties in these parameters, we show what 
might be expected for a survey with the same number den- 
sity as SDSS-II (n ~ 10" 4 /i 3 Mpc~ 3 ) but a volume twelve 
times larger, corresponding to the total Hubble volume (i.e, 
a sphere of radius c/Hq). 

To simulate the observer's point of view, we assume the 
£ v measurements to be the "data" points. We then find the 



best fitting models based on physical templates. Our first 
step is to perform an ideal test, where the template is the 
£ r mock mean signal. In other words, we are not concerning 
ourselves, at this point, with uncertainties beyond the AP 
effect. For example, this means we assume that we fully un- 
derstand the amplitude "bias", and dynamical ^-distortions. 
We consider this merely as a "proof of concept" of the anal- 
ysis, and in H4.3I take a more realistic approach by adding 
more unknowns. 

In i ]2.1l we describe the construction of the covariance 
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matrix, in which we take into account covariances between 
the statistics. When manipulating the template, the covari- 
ance m atrix is not varied, but rather fixed to the true cos- 
mology. ISamushia et al.l (|201 ll ) discuss the sensitivity of Cy 
to amplitude parameters and its insensitivity to shape pa- 
rameters. 

Our results are presented in Figures[3][S] Figure 0] shows 
the two-dimensional marginalized constraints on the Da — H 
plane. The contours shown correspond to 40, 68, 95, and 
99% confidence levels. The best fit values of each parameter 
and their respective la (68%) CL regions after marginalizing 
over the other are indicated by the purple crosses. We find 
that the best fit model recovers, to high accuracy, the true 
values of these parameters. The pink box corresponds to the 
fiducial cosmology. In the legend we display the calculated 
deviations from the true values in the 2D analysis (purple di- 
amond). We also note the marginalized uncertainty for (i.e, 
the analog of the diagonal elements in Fisher-Matrix anal- 
ysis). For the wedges we make use of first order correction 
terms C± 11(e), and in Appendix [Bl discuss their importance. 
For comparison, the wide orange contours display the con- 
straining power of the monopole on its own. The result is 
the well known D\/H degeneracy. 

The fitting range chosen here is 40 < s < 150 fc~ 1 Mpc 
in bin s of As = 5fc~ 1 Mpc. As discussed in ISanchez et al.l 
(|2008h and lShoii et all < |2009h ■ we find that the extra infor- 
mation included when analyzing larger ranges yields tighter 
constraints. 

We notice that in all cases the true parameters are re- 
covered to high accuracy, and the input incorrect cosmology 
is ruled out by over 3a solely by the AP effect on cluster- 
ing. That said, this level of precision is not expected when 
using this technique in near future galaxy samples, because 
marginalization over a larger parameter space would be re- 
quired. We discuss this more detail in H4.3I 

We do notice, though, differences in the results obtained 
by means of the two projection techniques. We emphasize 
that these are fair tests, as every step along the analysis is 
equivalent as much as possible. When analyzing real-space 
information (left plots), we see that Da is measured to 
similar accuracy with both techniques, with the clustering 
wedges yielding slightly smaller uncertainties. Interestingly, 
the differences in the recovery of H between the two meth- 
ods are larger where the clustering wedges yield a more ac- 
curate result, as well as smaller uncertainties. Throughout 
this study we compare these statistics, and find that the 
clustering wedges defined by A/j, = 1/2 perform better than 
£o,?2, for the most part, and motivate adding £4 to improv- 
ing constraints when using multipoles. 

In the plots on the right we compare the performance 
of these statistics in velocity-space. Most notable is the fact 
that the uncertainties in H increase substantially for both 
projection pairs compared to real-space results. Although 
the wedges yield tighter uncertainties in H, the multipoles 
generate slightly more accurate results. As for Da they both 
yield similar uncertainties, but clustering wedges yield more 
accurate results, and do not change much from real-space. 

When fitting for the A/j, — 1/2 clustering wedges, we 
take into account the intermixing terms C\\_± (Equations 1231 
and I24p . In Appendix [B] we explore their effectiveness and 
test results obtained by other wedge widths. We find that 
setting Cm 1 = (meaning radial wedge depends solely on 



H and transverse wedges solely on -Da), yields results that 
are less accurate and uncertainties that are underestimated. 
We also find that the C\^± terms are less important as one 
decreases A/i. 

In the top panels of Figure [S] we compare the results 
shown in Figure[3]with ones obtained when adding £4 to the 
multipole pair. The bottom plots corresponds to the same 
test but using a different choice for the fiducial cosmology 
(w = —0.9; notice the different location of the magenta box). 
The left plots are in real-space and the right in velocity- 
space. 

In all cases the true cosmology is recovered to high ac- 
curacy. As expected, the [^0,^2,^4] combination (dash blue 
lines) yields tighter constraints than when limiting the anal- 
ysis to [£0,^2] (solid black lines) or the Afi = 0.5 wedges 
(dot-dashed orange lines). Crosses indicate marginalized la 
results according to color, and symbols most likely 2D value. 
The real-space w = —0.9 result (bottom left plot) shows that 
the wedges do not generically yield tighter constraints than 
the monopole-quadrupole pair. In velocity-space, we notice 
that adding the £4 information in the multipoles yields a 
smaller correlation coefficient between H and Da- 

Another oddity is the fact that in real-space £4 is present 
at all, and can assist in improving constraints. We argue that 
it is probably an artifact of angular effects at large scales (see 
H4.1|) . We verify this by limiting our analysis to the range 
s = [20, 60] h~ Mpc, finding less of an improvement when 
using £4, as its angular effect is negligible at these scales. 
These suggest that including £4 is useful in breaking the 
H — Da degeneracy in velocity-space. 

To summarize the results of this section, we find that 
in the ideal case analysed here the true parameters are re- 
covered to high accuracy with both statistics. The input in- 
correct cosmology is ruled out by over 3a solely by the AP 
effect on clustering. We demonstrate here, for the first time, 
the power of the clustering wedges technique. We find that 
using Equations (|21|) and (|22[) to describe the AP distortions 
in the clustering wedges it is possible to recover the true pa- 
rameters as well as one can with the multipole expansion. 
This means that even a wide "radial" wedge is most sensi- 
tive to H while a "transverse" wedges is most sensitive to 
Da- We also demonstrate that adding £4 to the multipoles 
substantially improves constraints. 

4.3 Amplitude effects on uncertainties 

In the previous sections we simulate "ideal" observational 
tests, in which the correlation function template from which 
the models are constructed, is fully understood except for 
the AP effect. Here we address the fact that the observer 
does not have the luxury of knowing the correct signal a 
priori. 

To mention a few main concerns, one must understand 
the £ amplitude, shape, and various effects on the angular 
baryonic acoustic feature. The amplitude may be scale de- 
pendent for various reasons: non-linear tracer-matter bias, 
scale dependent dynamic disto rtion effects, ob server angu- 



lar effects, magnification bias (|Hui et al 



2007) and effects 



of non-gaussianities in initial conditions ( Dalai et al.ll2008T ). 

Here we take into account the effect of uncertainties 
in the amplitudes of the various statistics and leave shape 
effects for a future study. For this, we add amplitude param- 



© 0000 RAS, MNRAS 000, 000-000 



Improving measurements of H(z) and D\(z) by analyzing clustering anisotropies 11 




Figure 4. H — D\ joint constraints from AP shift tests when analyzing mock mean clustering ID statistics. Left panels: real-space. 
Right: velocity-space. The blue ellipses correspond to analysis of [§o,£2] (top panels) and A/i = 1/2 wedges (bottom). The orange bands 
are the results when analyzing only the go- The parameter space tested is $ =[D\,H], where amplitude remains fixed ("AP only"). 
The true cosmology of the simulations corresponds to [1,1], and the distortion cosmology tested when converting redshifts to comoving 
distances (pink box) corresponds to using = —1.1 instead of the true —1 value, and is the same for all tests. The even contour 
likelihood mock mean values are for 40, 68.3, 95.4, 99% CL. The most likely 2D values are purple diamonds, and the ID marginalized 
la results are the purple crosses. On the top of each plot we inscribe for each parameter <& the deviation from the true value in the ID 
marginalized results, as well as The constraints correspond roughly to a fictitious survey with a Hubble volume (R ~ c/Hq) and 

a galaxy density of n ~ 10~ 4 (/i/Mpc) 3 . 



eters to the H — Da parameter space. We assume that the 
amplitudes of the multipoles and clustering wedges are given 
by A s tat = -B s tatb 2 a|, where b is the (scale-independent) 
bias parameter of the tracer with respect to the matter, ag is 
the rms linear perturbation theory variance in spheres of ra- 
dius 8/i _1 Mpc, and _B s tat is the Kaiser squashing amplitude 
on large scales, which in real-space is unity. The subscript 
"stat" signifies the fact that we test for various statistics 
such as the wedges and multipoles. 

For the multipoles in velocity-space B 8 t a t is given by 



the usual iKaisej (1 19871 ) prefixes 

Bo = 1 + 2/3^ + 1/5/3 2 , (25) 

B 2 = 4/3/3 + 4/7/3 2 , (26) 

B 4 = 8/35/3 4 , (27) 

of the squashing parameter /3(z) = f/b ~ Q.m(z) ' 55 /b, 
where / = dln(Di)/dln(a), with D\ the standard linear 
growth factor and and a the scale factor. Technically, as the 
templates used here correspond to the T signal, when test- 
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Figure 5. Joint 2cr constraints on the H — D\ plane from AP shift tests when analyzing mock mean clustering ID statistics. Left panels: 
real-space. Right: velocity-space. The AP distortion tested is when using w different from the true value iifT = — 1 when converting 
redshifts to comoving distances. Top panels show results for w' D = —1.1, and bottom panels for w D = —0.9 (notice difference in input 
fiducial cosmology indicated by the magenta squares). The Parameter space analised is [Dp^,H\, where the amplitude remains fixed 
("AP only"). The dashed blue lines correspond to the constraints for the multipole combination [£o>?2,£4], solid black lines to [£o,?2], 
and dot-dashed orange lines to Afi = 1/2 wedges. Crosses indicate marginalized la results according to color, and symbols most likely 
2D value. The w° = —1.1 clustering wedges and [£0,^2] results are the same as the blue ellipses in Figure [4] Adding the hexadcapole 
information clearly improves constraints. These constraints correspond roughly to a fictitious survey with a Hubble volume (i?, ~ c/Hq) 
and a galaxy density of n ~ 10~ 4 (/i/Mpc) 3 . 



ing for the amplitude the models are based on the ratio of 
j4 sta t to its real value. 

When testing for bag we perform the analysis on wedges 
and multipoles. When testing for /?, however, this is a non- 
trivial task for the wedges. This is due to the volume aver- 
aged monopole term which makes it non-trivial to account 
for A s tat in configu ration space (se e the terms involving £ 
in equations 6-8 of iHamiltonl 1 199^ ) . In a k— space analysis 
this issue would be trivial, as it should be when building a 
generic model, which we do not do here. For this reason, re- 



sults shown marginalizing over (3 are limited to multipoles. 
The fiducial f} tested is calculated through the input f(z) 
and the linear b ~ 1.96 which is inferred from matching 
linear theory to the standard projected correlation function 
w p (s±) of the mocks. 

In Figure [6] we show joint constraints on H and Da us- 
ing the various statistics when marginalizing over fcerg (top) 
and [6<T8) 0\ (bottom). For the former case we assume perfect 
knowledge of the dynamic distortions, whereas for the lat- 
ter we add the uncertainty of the Kaiser effect (disregarding 
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Figure 6. Joint 2a constraints on H — Da from AP shift tests 
when analyzing mock mean velocity-space clustering ID statis- 
tics. The parameter space in the upper panel is <E> = [Da, H, bag]. 
and <!> = [ Da , H, bag , jS\ in the bottom. The analysis of the re- 
sults shown here is similar to that performed when producing 
the results shown in Figure [5] only marginalizing over amplitude 
parameters. H — Da joint constraints for the multipole combi- 
nation [£o, £2, £4] are shown in dashed blue lines, [^0,^2] are in 
solid black lines, and Afi = 1/2 wedges are dot-dashed orange 
lines. Crosses indicate marginalized la results according to color, 
and symbols most likely 2D value. The AP distortion applied 
here is using u> FID = —1.1 (orange box) instead of the true value 
— 1 when converting z to comoving distances. These constraints 
roughly correspond to a fictitious survey with a Hubble volume 
(R ~ c/ Hq) and a galaxy density of n ~ 10~ 4 (/i/Mpc) 3 . 



velocity dispersion). The obtained constraints can be com- 
pared with those of Figure [S] and show a clear degradation 
of the uncertainties on H nearly by a factor of two and a 
slight degradation of the uncertainties on Da with respect 
to the results obtained without uncertainties in the ampli- 
tudes. In i )4.4l we study the effect of varying the range of 
scales included in the analysis. These plots show fair com- 
parisons between [£0,^2], [£0,^2, £4] and A/x = 0.5 clustering 
wedges. 

We see that the true cosmology is recovered to high 
accuracy, but our choice of incorrect cosmology w = —1.1 is 
ruled out by [£0,^2] at only 2a (top) and ~ 1.5a (bottom), 
due to the marginalization over the amplitude parameters. 
The importance of the first finding is that multipole AP 
Equations (|19p and (|20p . hold to a very good degree. We 
obtain similar conclusions for Equations (|21|l . (|22|l and (|23|l 
regarding the clustering wedges. As in the test performed 
we find an increase in the uncertainties of H at the 



10% level when going from real- to velocity-space. 

To conclude these tests we find that increasing the al- 
lowed range of the amplitude parameters degrades the con- 
straints, but retains accuracy in recovering the true val- 
ues of H and Da- When comparing between the different 
statistics our conclusions of this test are similar to those ob- 
tained in i]4.2l We demonstrate that £4 improves constraints 
substantially compared to those obtained when limited to 
monopole-quadrupole. 

4.4 Wedges or multipoles? 



In £14.21 and £|4.3l we analyse a particular case where the con- 
straining power on H of the A[i = 0.5 clustering wedges 
outperforms the monopole-quadrupole pair, and see sub- 
stantial improvement when including the hexadecapole (see 
Figure |3J. However, in these analyses we limit the range 
of scales to 40 < s < 150 h~ Mpc. Here we generalize 
these tests by varying the minimum scale included in the 
analysis s m i n while keeping the maximum scale fixed to 
Smax = 150 /i _1 Mpc, and find interesting trends for the con- 
straining power of the various statistics. All results given 
here are in velocity-space. 

Our findings are summarized in Figure [7] in which we 
show the uncertainties on H and Da as a function of s m i n . 

Results for an " AP effect only" are displayed in Figure 
[7^, which generalize the top right plot in Figure [S] " AP 
only" refers to tests in which we fix amplitude parameters 
to their true values and test only for H and Da- When 
marginalizing over bag we obtain results displayed in Figure 
[7p, which generalizes the top plot of Figure [6] Figure [Jj; 
shows the results obtained when marginalizing over both 
bag and /3, and is a generalization of the bottom plot of 
Figure [6] No priors are assumed for bag or j3. 

In all cases we find that adding information from £4 to 
the multipole analysis improves the obtained constraints on 
H and Da by a substantial amount. 

The improvement in the constraints as smaller scales are 
included in the analysis emphasizes the fact that although 
the baryonic acoustic feature is essential to perform the AP 
test, one can extr act more information by analyzing the full 
broad shape of £ (|Sanchez et al.ll200Sl ; Ishoii et al.ll2009l ). 

In all parameter spaces tested, the slope of anismin) 
(thick red) is steeper than that of a d a (thin blue) , in- 
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Figure 7. Relative uncertainties ajf/H and o"d a / d a as a function of the range analysed [s m i n , 150 h~ 1 Mpc] (given in [%]). The 
different symbols, indicated in the legend, represent the various ID statistics used. Thick red symbols are for an, and thin blue for <j n,. 
(a): " AP only"- [H, Da] with fixed amplitude, (b): [H, Da] when marginalizing over berg without prior, (c): [H, Da] when marginalizing 
over bag, /3 without prior. Slopes indicate the improved constraints obtained by using information from the broad band of £. The dotted 
lines indicate an improvement of 0.8 from <r 4 ,/<I>(s max = 80 h~ 1 Mpc) (when fitting for [£o.£2)£4]i ^ = [H, Da]), and the dashed lines an 
improvement of 0.5. The AP distortion applied here is using ui FID = —1.1 instead of the true value —1 when converting z to comoving 
distances. 



dicating that the broad shape of these statistics is more 
sensitive to H. The fact that the "AP only" test (Fig- 
ure [7^,) has steeper slopes than the others is expected. In 
this case each data poin t posse s a high constraining power. 
IPadmanabhan fc White! (|2008l ) explain that marginalizing 
over the amplitude is somewhat degenerate with the AP ef- 
fect and hence leads to a reduction of the slope when adding 
bag, (Figure (TJd) and even more when adding /3 (Figure [7t). 
Nevertheless, even when marginalizing over bag, we notice 
an improvement of factor two in the uncertainty on H and 
a factor of ~ 1.25 in a d a when using the full information 
up to 40 /i _1 Mpc, with respect to setting s m i n = 80 /i~ 1 Mpc 
(i.e, focusing on the scales around the baryonic acoustic fea- 
ture). 

The answer to the question of which statistics should be 
preferred is not simple as it appears to depend on the range 
of scales used. For high values of s m - m the results obtained by 
means of the clustering wedges (diamonds) outperform those 
obtained with the [£0,62] combination (triangles). However, 
we notice the different slopes of o"D A (smm) and a\Tf(s m i n ) 
such that the multipole pair should be preferred when a 
large range of scales is considered (i.e. low s m i n ). 

A puzzling result is the fact that a d a (s min) and 
o\ff(smin) have a steeper slope for the case of the [£0,62, £4] 
combination than for the clustering wedges. The (J_ff(s m i n ) 
results indicate that the A/x = 0.5 wedges are preferred at 
Smin > 50 ft - Mpc (although [f 0,^2, £4] is preferred to de- 
termine -Da). This might be explained by angular effects 
causing higher multipoles, which we have not corrected for 



properly with the multipoles, but happen not to affect the 
wedges, which are combinations of all multipoles. 

For this analysis we have chosen w fid = —1.1 of our 
fiducial cosmology (instead of w tTUC — —1). When testing 
for w Rd = -0.9 we find similar trends but with varying 
crossover points. For example, in the "AP only" case, the 
clustering wedges are preferred over the [£o,f2,f4] combi- 
nation in determining H at s m i n > 55 /i -1 Mpc instead of 
> 45/2T 1 Mpc. 



5 PREDICTED CONSTRAINTS ON H AND 
D A FROM BOSS 

Here we apply both the multipole and wedge techniques to 
the 32 realizations of the Horizon Run mock galaxy cata- 
logues described in H2.2I which serve as simulated realiza- 
tions of the final BOSS volume. We investigate the accuracy 
with which H and Da can be obtained for a survey covering 
10,000 square degrees in the redshift range 0.16 < z < 0.6 
with density n ~ 3 ■ 10~ 4 /i 3 Mpc -3 . For full details of the 
mocks please refer to £12.21 Analyzing 32 realizations, we ob- 
tain results for the full sample, and of for two subsamples 
split at z = 0.45. This redshift split is motivated by the 
so-called targeted LO WZ and CMASS BOSS galaxies (see 
lEisenstein et alj|201ll and Padmanabhan et al. (in prep.)). 

As in previous sections, when building our templates, 
we gradually increase our parameter space, starting from 
solely the AP effect and adding amplitude parameters. For 
each case we analyse both real- and velocity-space and ob- 
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Figure 8. BOSS mock galaxy mean (0.16 < z < 0.6) H — Da joint 2a constraints from AP shift tests, when analyzing velocity-space 
clustering ID statistics. The parameter space in the left panel is [ D\,H ,bag] . and [ D\,H,ba g,/3] in the right. H— D\ joint constraints for 
the multipole combination [foj?2;£4] are shown in dashed blue lines, [£cbC2] are in solid black lines, and Afi = 1/2 wedges are dot-dashed 
orange lines. Crosses indicate marginalized la results according to color, and symbols most likely 2D value. The AP distortion applied 
here is using u> FID = —1.1 (orange box) instead of the true value —1, when converting z to comoving distances. These constraints 
correspond to an analysis of the broadband s = [40, 150] h~ 1 Wlpc, without priors, and assumes all shape effects are fully understood. In 
both [5oi^2,?4] cases the true cosmology is recovered with excellent accuracy, and the incorrect input H, Da is rejected by ~ 2a. We 
clearly sec a degradation in constraining H when including marginalization of f3. 




Figure 9. BOSS mock galaxy mean (0.16 < z < 0.6) joint and marginalized constraints from AP shift tests, when analyzing velocity- 
space clustering [£oj£2,£4]. (These results are the same as those displayed in Figure \E\ meaning analyzing region s = [40, 150] ft Mpc). 
The parameter space in the left panel is [ Da,H ,ba$] . and [ Da,H ,bag,p] in the right. The joint constraints are for CL: [40, 68, 95, 99]%. On 
the top of each ID likelihood plot we inscribe for each parameter $ the marginalized la result, also summarized in Table IT1 The orange 
boxes demonstrate the fiducial parameters input into the analysis (in the case of H and D\ this is the AP effect using u> FID = —1.1 
instead of the true —1 value when converting redshifts to comoving distances). In all cases the true cosmology is recovered with excellent 
accuracy, and the incorrect input H, Da is rejected by ~ 2a. We clearly see a ~ 30% degradation in constraining H when including 
marginalization of ft, and a ~ 15% degradation in Da- 
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tain the corresponding results as expected from one BOSS 
volume. 

On a technical note, the mock "data" and templates 
are based on the mock mean of 32 Horizon Run mocks. To 
obtain a stable, invertible dj we use the 160 LasDamas 
mocks, normalized by the variance of the 32 Horizon Run 
mocks. 

Our 0.16 < z < 0.6 results are shown in Figures [8] and [9] 
for a particular range 40 < s < 150 h~ Mpc. The left plots in 
each figure are results for the parameter space [H, Da, bag] 
and the right plots for [H, Da, bag, 0\. Figure [8] compares 
the statistic combinations in the H — Da plane. In Figure [9] 
we display one and two-dimensional likelihood functions for 
the above three and four parameter spaces when using the 
[£o, £2,^4] combination. Marginalized la uncertainty values 
are indicated on the ID likelihood panels. In Table [T] we 
summarize the la uncertainty values obtained for the vari- 
ous statistic combinations, parameter spaces, choices of s m i n , 
and sample analysed. 

In all cases investigated we find that the true cosmol- 
ogy (H, Da) is recovered to high accuracy (much better than 
1%). We report that the velocity-space results yield a notice- 
able increase in uncertainty of H compared to real-space re- 
sults. Below we will limit our explanations to velocity-space 
results. 

When allowing the amplitude parameters to vary, we 
find that uncertainties in the H — Da plane are degraded, 
as expected. When adding /? in velocity-space, we notice an 
increase in the uncertainty of H by a fraction range of 1.1 
to 1.6, depending on the scales included in the analysis. For 
Da the fraction changes by 1.01 to 1.10. 

As in the previous sections, we find in all cases that 
adding the £4 information imporves the constraints obtained 
by means of the multipoles. We see an improvement of a 
factor 1.2 to 1.35 in H and between 1.15 to 1.25 i n Da - 
This is comparable to results found in iTaruva et all i|201ll) . 
who perform a Fisher Matrix analysis. 

As for the clustering wedges, we find that they out- 
perform the monopole-quadrupole pair in H while giving 
similar constraints on Da- Interestingly, in this case the 
A/i = 0.5 clustering wedges measure H with a similar accu- 
racy to that of [£o, £2, £4]. 

As to values expected from BOSS, we find that, assum- 
ing that P is fixed and the broad shape understood down 
to 40ft _1 Mpc, our best 0.16 < z < 0.6 constraints obtained 
are: AH/H ~ 1.52% (similar for wedges and [£0,^2, £4]) and 
ADa/Da-1% ([£0,6, £4]). 

When splitting to the two subsamples, using 
[Co, 6, &] at 40 < s < 150/i -1 Mpc, NEAR yields 
[AH/H, A D A / D A }^ [2.13,1.46]% and FAR [1.85,1.41%]. 

ISchlegel et all i|2009l ) use Fisher-Matrix analysis to ob- 
tain much more optimistic estimates, than those shown here, 
even though they focus on the baryonic acoustic wiggles 
in P(k), and not the broad shape and do not use the £4. 
This is probably due to the fact that they assume the "re- 
construction" of the feature which, if applicable without 
introducing bias, shou ld improve the obtained constraints 
l|Eisenstein et al]|2007l ). We have not applied the technique 
on the mocks. 



6 DISCUSSION 

The purpose of this study is to investigate possible ways to 
break the H — Da degeneracy by including information in 
the anisotropy in the £ (fj,, s) plane produced by geometrical 
redshift distortions. 



6.1 Relating our analysis to previous studies 

This concept has bee n studied in the f ull 2D P(p,, k) plane by 
iHu fc Haimanl (|2003h . I Wagner et all (|2008l ), and lShoii et"atl 
1 20091 ). who investigated the power of using the baryonic 
acoustic feature to determine the equation of state of dark 
energy. In practical terms, however, there are a few difficul- 
ties in applying this approach on real data, namely the low 
S/N of the measurements in the full fi—k plane, the practical 
problems related to estimating accurate covariance matrices 
for them, and the difficulties in constructing realistic models 
that take non-lin earities into account. 

Following iPadmanabhan fc White! (|2008ft and 
ITaruva et alT (|201Ch . we break the H — Da degener- 
acy by using projections of the fj, — s plane, which have the 
advantage of a higher S/N, while preserving much of the 
essential information. As near future surveys will provide 
fairly noisy P{n, k) planes, we find the projection approach 
more useful in the short term. 

These last two studies focused on the monopole- 
quadrupole (or "multipole") pair in fc— space. We demon- 
strate similar results for the first time in configuration space, 
and introduce an alternative method in the form of cluster- 
ing wedges £(A/z, s), and compare their constraining power 
to that of the multipoles. 

The projection approach also simplifies covariance is- 
sues, as one resorts to a much smaller covariance matrix, 
which is more likely to be invertible and stable, when using 
a reasonable number of mock real i zation s. An alternative 
method suggested by ITaruva et ail (|2010h is to ignore non- 
Gaussianties by using a linear Cov^/(fc) based on an ana- 
lytic model. While this might be a fine approach for simple 
estimates, when analyzing real data one should take into ac- 
count observational effects, most straightforwardly achieved 
by mock realizations with a similar window function, such 
as the mock catalogues produced by the LasDamas group 
(McBride et al.; in prep), and used here. 

In th is study we analyse geom etric effects (or AP ef- 
fects after lAlcock fc Paczvnskl 19791 ) in clustering. We study 
the effect of using an incorrect value for the dark-energy 
equation of state w, when converting redshifts to comov- 
ing distances, causing slight shifts in the inferred positions 
of the galaxies in respect to the real positions. Instead of 
the true w = —1 value we use —0.9 and —1.1, which is 
within the allowed re gion for this paramet e r according to 
current observations (IKomatsu et al. 1 120091 , ISanchez et al.l 
20091. iPercivaT et al. I l2010l . iReid et al.ll2010t ). It is interesting 
to note that the —0.9 shift (which causes larger dilations and 
warps) yields slightly, but noticeable, larger uncertainties in 
both real- and velocity-space (to see this compare the cor- 
responding results in Figure [5| . It would be interesting to 
see if this is a systematic trend increasing with dilation and 
warping, or if we obtained this by chance. Although we ob- 
tain similar absolute results in H and Da with the different 
AP effects, this does point out to a possible systematic in 
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Table 1. BOSS predictions based on velocity-space Horizon Run mock galaxies" 
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a We use the dj based on 160 LasDamas (SDSS-II LRG) mocks normalized by the Ca of 32 Horizon Run (BOSS) mocks. 
b All BOSS mocks arc volume limited and cover 1/4 of the sky with n ~ 3 ■ 10~ 4 fe 3 Mpc~ 3 . For 0.16 < z < 0.6: V = 4h" 3 Gpc 3 , 

0.16 < z < 0.45: 1.8 /i" 3 Gpc 3 , 0.45 < z < 0.6: 2.2 h~ 3 Gpc 3 . 
c la meaning 68.4% CL when all other parameters marginalized over. 
d "AP only" means amplitude is fixed and we test for H, Dj±. 
e Analyses of £4 assume it can be measured and modeled for. In BOSS, we do, however, expect low S/N at baryonic acoustic 

feature scales. 



the estimated ucertainties. This could be tested by applying 
different fiducial cosmologies on actual data, and comparing 
final results. 

iPadmanabhan fc White! (|200ct ) examine much larger 
warps (e > 1 compared to our more realistic |e| ~0.003) 
and show a trend of increasing uncertainty cr e with increas- 
ing e for e > 2 (see their Table 1). Their argument for using 
large e is that warping is degenerate between dynamic and 
geometric z-distortions, w h ere th e former clearly dominates 
the latter. iBallinger et al.l (|l996l ) suggest that the degener- 
acy between dynamic and geometric distortions may be re- 
solved by using measurements at various redshifts, as they 
are affected differently. 

W e perform similar analyses as IPadmanabhan fc White! 
(|200ct ) on the multipoles up to a few technical differences, 
which should not affect the results. The first difference is 
that they warp the box of their simulation, that is, distort 
the positions of the particles according to given values of 
a and e, while we imitate the observer's point of view by 
assuming an incorrect cosmology when converting redshifts 
to comoving distances. Second, we constrain both dilation 
a and warping e parameters simultaneously, while they as- 
sume q is constrained by the monopole independently, and 



this information is then combined with the quadrupole to 
constrain e. 



6.2 Modelling issues 

In this study we avoid modelling issues, for the most part, 
by using the true mock-mean signal as a template. By do- 
ing so we assume all parameters and effects are known 
except for the AP effect (H, Da), and test the effects of 
marginalizing over amplitude parameters bog, /3. As this as- 
sumes ideal conditions, and we do not test shape effects 
and non-linearities, we consider the tests performed here 
merely as proofs of concept. This means that we show that 
Equations (1171) (1201), init i ally in troduced in fc— space form by 
IPadmanabhan fc White! (|200cf ). and our equivalent version 
of the clustering wedges, Equations (121 1) (|23|) . accurately de- 
scribe the AP effect on projections of £(fi, s), and can be 
used to obtain constraints on the values of H and D\. 

In our tests on mock catalogues we show that the AP 
effect does not introduce substantial amplitude or shape ef- 
fects. Minute amplitude deviations are seen in the erg /aJ RlJE 
column in the bottom left plot of Figure |§J as deviations of 
the results from unity. This amplitude test shows that the 
AP effect does not introduce an amplitude bias. 
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We find that when marginalizing over /3 without priors, 
H constraints are substantially degraded. This is clearly seen 
in Figures [Sj [9] (degradation by factor of ~ 1.5) and Tabled] 
(For a generalization of this effect as a function of range of 
analysis see Figure]?]) The uncertainties on j3 are fairly large, 
too . The latter c ould be decreased by performing in parallel 
thelKaiscr ( 1987 ) quadrupole test based on the squashing ef- 
fect (see iTocchini- Valentini et al]|201ll for a newly proposed 
method to reduce uncertainties on /3 through the quadrupole 
test). 

In order to use the AP correction in practice on the 
broad shape of anisotropic clustering projection, or even if 
focusing only on the baryonic acoustic feature, a realistic 
model based on physical principles should be used as tem- 
plates for these statistics. Of special concern is understand- 
ing the distortions of the baryonic acoustic feature itself. 
For example, comparing the linear theory for £2 of Figure 
Cll with the results from the mock catalogues of Figure [3] 
we notice that the baryonic acoustic feature is distorted. In 
this case it appears to change from a dip to a bump. 

The only studies that the auth ors are aware o f that 
atte mpt to resol v e th is issue are iTaruva et all (|2009l ) 
and ITaruva et al.l (|2010T >. Following a model that includes 
veloci ty-dispersion decompression given in IScoccimarrol 
they improve the standard perturbation theory 
velocity-space power spectrum. They show significant im- 
provement compared to linear theory and previous attempts 
(see references within). They conclude that density and ve- 
locity terms need to be improved, as well as scale-dependen t 
and stochastic effects of galaxy bias. ISamushia et all (|201lT ) 
also show fairly good fits to the LasDamas mock catalogue 
using a more simplistic approach (see their Figure 11). 



6.3 More practicalities 

An approach that could potentially reduce non-linear ef- 
fects is the reconstruction techniqu e of the baryonic acous- 
tic feature. lEisenstein et al.l ([2007) proposed to reconstruct 
the monopole to its original linear form by reversing the 
displ acements of galax i es using t h e Zel 'dovic approxima- 
tion l |Zel'Dovichlll97Ch . ISeo et al.l (|201Ch demonstrate for 
dark matter particles in real- and velocity-space that the 
non-linearities can be corrected for to a high degree, and 
iMehta et al.l (|201lh have recently reported that reconstruc- 
tion should work for biased matter tracers in an unbiased 
manner. Focu sing on the feat ure, and ignoring effects of 
the full shape, ISeo et all l|201(t ) show that the a parameter 
(~ D\/H) can be reproduced accurately in unbiased fash- 
ion. The fact that their velocity-space results at low redshifts 
show a larger s catter than in real-space (see Figures 3-5 in 
ISeo et all I2010T 1 indicates that there is still information in 
the higher multipoles, especially the quadrupole. Nonethe- 
less, this is encouraging in the matter case, and it would be 
interesting to see if the remaining quadrupole after recon- 
struction would be useful to constrain cosmology. If the dy- 
namic £2 could be eliminated, however, Equation (|20[) yields 
the simple relation t$ = 2ed£o/dlns, and hence one needs 
to model only for the monopole, which in this case would be 
(very close) to its original linear form. 

One needs also to consider the fact that dynamical 
distor tions are degenerate with geometrical. iBallinger et al.l 
(|l99rj ) argue that, although breaking this degneracy might 



be impossible at one given redshift, analyzing various epochs 
could break degeneracies as both effects evolve differently. It 
would be interesting to see if this could be done in practice, 
where in reality one might be using different sets of tracers , 
as well as amplitude bias evolution of the tracers (|Frvlll996h . 



7 CONCLUSIONS 

We demonstrate that by correcting for the geometric effects 
of ID projections of the clustering s), it is possible to 
constrain H, Da to high accuracy We perform tests on the 
commonly used monopole-quadrupole pair as well as an al- 
ternative basis in the form of clustering wedges, introduced 
here for the first time. By doing so, we pro ve that the ge- 
ometrical effects (|Alcock fc Paczvnskilll979T ) are accurately 
described by Equations (|21|) - (|23| l. which illustrate that even 
a wide "radial" wedge /i > 0.5 is mostly sensitive to H and 
a "transverse" wedge (/j, < 0.5) is sensitive to Da, up to 
small intermixing corrections terms Cm j_. 

Throughout this study we use both analytic formu- 
lae and realistic mock galaxy catalogues to compare the 
constraining power of the wide A/j, — 0.5 wedges with 
the previously pr oposed multipole statistics, the monopole- 
quadrupole pair C lPadmanabhan fe White] I2008I ). Our main 
findings are: 



(i) The lAlcock fc Paczvnskl (| 19791 ') effect in £ is very well 
described by Equations (0 and pOjl. (|2T ]l -p3 ]l (e.g, see 
Figures H EH. 

(ii) Adding the hexadecapole in the multipole analysis 
improves constraints on H by a factor of 1.2 — 1.35 and 
improves Da measurements by 1.1 5— 1.25 (Figures [5H7}. In 
a recent studv lTaruva et al.l (|201ll ) have used Fisher matrix 
analysis and obtained similar results. 

(iii) The clustering wedges serve as an alternative basis 
(see Figure ID1|) containing much of the same H, Da in- 
formation as the standard multipole projection (e.g, Figure 
0). 

(iv) Limiting the analysis to the baryonic acoustic feature 
region, the Afj, = 0.5 clustering wedges can slightly out- 
perform the [£0,^2, £4] combination in constraining H (see 
Figure [7}. This might be due to the fact that higher multi- 
poles need to be taken into account. See t|4.4l for a detailed 
discussion. 

(v) Constraints on H and Da can be substantially im- 
proved by analyzing the broad shape of £ (see Figure [7] and 
discussion in £14.41). 



The improved constraining power of clustering wedges 
might be explained by the fact that they contain informa- 
tion from higher order multipoles (see Equation [S}. We also 
argue that at "low" redshift (e.g, (z}=0.33 as tested here) 
angular effects introduce even higher order multipole con- 
tributions which might have an effect at baryonic acoustic 
feature scales, when correcting for the AP effect (see Figure 
1101 and Appendix ID] for details). 

To predict how well the ongoing BOSS survey will do, 
we analyse mock galaxies of the volume and number density 
expected for this sample. Although we predict improvements 
in constraining H by using information on £4, this assumes 
both measuring a signal and being able to model it. In prac- 
tice, BOSS should yield low S/N at baryonic acoustic fea- 
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Figure 10. Mean high order multipolcs from the set of 160 Las- 
Damas m ock catalo g ues co mpared with linear theory pred i ctions 
based on iHamiltonl dl992T) (dashed line) and iMatsubaral II2004T ) 
(solid lines). The variances correspond to the uncertainty in the 
mean measure ments. Large ang le effects (which are included in 
the formulae of Matsubara 2004) reduce the amplitude of £4 and 
cause the apparent £6, £s results. 



ture scales, meaning that the use of £4 would only be possible 
at smaller scales, where velocity-dispersion effects which are 
difficult to model are more important. 

Because we test for only a few types of ID projections, 
there is room for further investigation for an optimal statis- 
tic for extracting cosmological information from clustering 
through the AP test. Nonetheless, having various types of 
statistics yielding consistent results provides an important 
tool to test for systematics. 

We suggest that a Fisher Matrix analysis could serve 
as an analytic method to test for an optimal method, to be 
followed up with mock catalogue tests. 

Once the above are achieved, Figure [2] demonstrates 
that degeneracies between Qmo, and the dark energy 
equation of state w within the H— Da plane may be resolved 
by using clustering at high redshifts (z > 2; e.g, through 
Lyman-a forest, and 21 cm measurements) in a complemen- 
tary fashion to other observations (e.g, the temperature fluc- 
tuations of the CMB). At lower redshifts the degeneracies 
between these parameters is quite large at a single redshift, 
but might be resolved by measuring H(z) at various red- 
shifts. 

In the coming years a variety of new large volume 
galaxy surveys will measure the large-scale clustering pat- 
tern of the Universe with unprecedented precision. Inves- 
tigations of techniques su ch as the reconstruct i on of the 
baryonic acoustic featu re (|Eisenstein et all 120071 ; ISeo et all 
l201fj ; [Mehta et al"1l201ll . Padmanabhan et al (in prep.)) sug- 
gest possible substantial improvements on constraining the 
cosmic evolution out of the information from these surveys 
and hence on our understanding of dark energy. That said, 



these measurements will always be bound to the AP effect. 
The ID clustering projection techniques discussed here will 
be essential to ensure that the full potential of the informa- 
tion contained in the AP effect will be realised. 
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APPENDIX A: HEXADECAPOLE TERMS 



In H3.4l we give analytic expressions for the AP effect on the 
monopole and quadrupole. Here we extend this treatment to 
take into account the hexadecapole £4 contribution. In this 
case Equation Q20p becomes: 



£3 to = £a(aa) + e 2-^-( + -&(as) + - 



4 <%2(S 



dln(s) 



7dln(s) 



20,. , 4 d£ 4 (s) 
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In our analyses of the AP effect on the £2 measurement 
from the mock catalogues, we find the £4 and d^/dlns cor- 
rections negligible. Limiting the hexadecapole to I <= 4 
contributions we obtain: 

£4 to = £ 4 (as) + e — , , - =r£ 4 (as) + 



35 dln(s) 77 



154 dln(s) 



(A2) 



where a higher £ would be required for completeness. 



APPENDIX B: TESTING WEDGE WIDTHS 
AND INTERMIXING TERMS 

In Figure IBTI we explore the effectiveness of the wedge cor- 
rection terms Cm 1 (Equation [23} • These results should be 



compared to those shown in the bottom plots of Figure |4j 
The short-dashed blue lines are the 2a results shown in Fig- 
ure[4] meaning when including C\\.± terms. The thick single- 
dot-dashed black lines are the 2a results from the same tests 
where we set C^x = 0. Two interesting differences are ap- 
parent. The most obvious one is the small bias relative to the 
true cosmology which produces a shift in the contour lines 
(although within the la region). Interestingly, our second 
observation is that it yields apparently tighter constraints 
than the corrected method (colored contour). These obser- 
vations are noticeable in both spaces. 

In Figure IBTI we also investigate other choices of radial 
and transverse wedges and give 2a results for: Afi = 0.125 
and 0.25, where we do not use correction terms (i.e, Cm ±=0). 

Our results clearly show, as expected, that using more 
information, i.e., through wider wedges yields tighter con- 
straints. We also notice that decreasing A/j, reduces the bias 
in the obtained constraints on Da both in real- and velocity- 
space. In velocity-space the bias is improved also in the H 
direction. We have not investigated the constraining power 
of multiple thin wedges, but rather that of the two extermes 
(i > Afi and /i < Ajx. 

In both spaces we notice that the narrower clustering 
wedges yield contours with reduced correlation coefficients 
between parameters H and Da- This is expected due to 
weaker intermixing terms with decreasing A/i. This is most 
prominent in velocity-space where the contours sharpen to- 
wards the line-of-sight. In real-space the ellipticity appears 
to decreas with A^i. 



APPENDIX C: PROOF OF CONCEPT: USING 
ANALYTIC PREDICTIONS 

Here we demonstrate the applicability of Equations (|19|l 
((23} to reproduce the distortions in the various multipoles 
and clustering wedges introduced by the use of an incorrect 
cosmology when transforming the observ ed redshifts todis- 
tances. For this we use the formulae of iMatsubaral (|2004T l 
for the velocity-space £(/i, s). We use their equation 1 which 
yields the two-point correlation function £(zi, z%, 7) for pairs 
of objects located at redshifts z\ and Z2 separated by ob- 
server angle 7. Applying Equation (|3]), and basic Euclidean 
comoving geometry, we project this result around Z\ to ob- 
tain a s) plane, where the line-of-sight direction is de- 
fined as the bisecting vector from the observer to s (i.e., the 
same as we do with the mock catalogues). We then project 
£(/z, s) into the various multipole and wedge components. 

For the true signal ( r ( we use the true input P(k) 
cosmology in Equation @. For the distorted signal ^ we 
apply the AP effect in which we apply the same exact pro- 
cedure, except that in Equation ((3} we assume an equation 
of state w = —1.1 instead of the true —1 value. For the 
shifted signal ^ we apply the scaling given by Equations 
(|19|l and (|20|l on the £ r signal of the monopole, quadrupole 
pair (and also on £4), and Equations (|2ip. ()22p . and (|23|l for 
the A/i = 0.5 wedges. 

Our results are shown in the upper panels of Figure 
ICll for the multipoles (left) and clustering wedges (panels), 
rescaled by s 2 for clarity. The bottom panels compare the 
predictions of Equations (fl9 )l -p3 )l with the £ D (or £ AP ) 
results. A perfect description of the AP effect would result 
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Figure CI. Analytical velocity-space clustering projections (based on lMat subara 2 004h with and without the AP effect. Left: multipolcs 
(corresponding to Equations 1191 l20ll . Right: A/i =1/2 wedges (Equations I2H231) , In the top panel of each plot are the ID projections, 
and in the bottom panels the difference of each result with the DISTORTED signal (labeled as AP) normalized by the uncertainty of 
one realization. As indicated in legend, solid lines are the true projection signals (T), the long-dashed lines are the AP signals (T>). We 
apply the AP shift on true signals to obtain the triple-dot-dashed lines ( S). A perfect shift would yield a null result in the bottom panels. 
The symbols represent applying only zeroth order corrections. For multipoles this means applying only isotropic dilation correction (a) 
and no warping terms (e). We see that the warping is crucial for correcting £2 but negligible (although visible here) for the monopolc. 
For the wedges zeroth order means applying only H ratio corrections for radial wedge, and applying Da ratio for transverse (meaning 
no C|| x corrections). The AP distortion applied here is using «) FID = —1.1 instead of the true value —1 when converting 2 to comoving 
distances. 



in zero values along the solid line. By comparing results of 
(&tat ~ 6rtat)/o*Btat and (£j at - f£at)M s tat we verify that 
Equations (|19[) (|23l) accuratedly describe the AP effect and 
can be used to correct for it. Here "stat" means the various 
ID statistics investigated: £0,2,1 

We also examine the importance of the first order e 
correction in each projection relative to the zeroth order. In 
the multipoles the zeroth order term is denned as applying 
only the isotropic dilation parameter a, meaning assuming 
tu( s ) = (<(as). In the wedges the zeroth order means apply- 
ing only the H ratio in the radial, and Da in the transverse. 
This means assuming Cm x = in Equations (|21[) and (|22[) 
for A/i = 1/2. 

The first order correction in all projections is denned 
here by adding the warping term e which include derivatives 
of the projections. 

The improvement in the monopole is subtle, and negli- 
gible for any practical case. In £2 the e term is, of course, es- 
sential for describing the AP shift. In the A/j, — 0.5 wedges, 
interestingly, the first order corrections yield slightly bet- 
ter results, meaning that the radial wedge is most sensitive 
to H and the transverse on Da- In Appendix [Bl we show, 
using mock catalogues, that the C(e) intermixing terms are 
essential to use in the wedge AP correction technique to 
yield unbiassed H, Da results, and that without them the 
uncertainties in these parameters are underestimated. 



APPENDIX D: PROJECTIONS IN PRACTICE: 
RELATING WEDGES TO MULTIPOLES 

In Figure ID II we put Equation (|lip to the test both ana- 
lytically and using mock galaxy catalogues by comparing 
A/j, — 0.5 wedges which are calculated directly to those ap- 
proximated by the multipoles. 

All uncertainty bars in the mocks are for the mock 
mean. All results (analytical and mocks) assume measure- 
ments of galaxies at (z) ~ 0.33 with a bias of b ~ 2 with the 
same flat ACDM cosmology. 



Dl Relating multipoles and wedges: analytic 
treatment 

In this section we test the relation between the clustering 
wedges and the multipoles by using theoretical p redictions 
for £( /x, s). For this we use the pres c riptio ns of iHamiltonl 
|l992T ) (Figure |D"Th) and IMatsubaral (|2004l ) (Figure iDlb). 
Both of these formalisms take into account only the Kaiser 
squashing effect in the plane-parallel approximation (and 
no velocity-dispersion), but the latter also includes wide- 
angle effects as well as linear growth evolution between 
points «i and Z2 (see Appendix [C] for a brief explan ation 
on manipulat i ons of Equation 1 of IMatsubaral [20041 ) . For 
the IHamiltonl | 19921 ) prescription we calculate £0, £2, and 
£4 directly by using their Equations 6-8, where the matter 
correlation functio n is calculated from th e output of CAMB 
jLewis et al.ll2002l ). For the formulae of IMatsubaral (|2004T l 
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Figure Dl. Panels (a),(b),(d) show the difference between clustering wedges inferred from multipolcs £^,cdgc X ' ant ^ those 

calculated directly, £™d C c' normalized by the uncertainty of one mock realization. Cljfj L. indicates the multipolcs [£] used in the ap- 
proximation. Panel (a) show results for analytic formulae in iHamiltonl j 19921 ), and (b) Klatsubaral d2004T) . respectively. Including the 
contribution from higher order multipolcs improves the results with respect to those obtained when only the monopolc-quadrupolc pair 
is used. Panels (c),(d) show similar results for the LasDamas mock catalogues, where (c) shows the actual statistics. At large scales, 
£o and £2 are sufficient to accurately describe clustering wedges obtained from the mock catalogues, and higher order multipoles are 
required for s < 30/i _1 Mpc. 



we use the same matter £(r) to compute f;(^i,s)n and in- 
tegrate to obtain multipoles using Equation ©Ij In each 
case, the "directly" calculated wedges are computed by in- 

4 In practice, £([i,s) requires very fine /i bins (actually, in 6). 
Here we show results for 89 ~ 5 ■ 10 — where for e.g. 68 = 1°, 
there are binning effects. 

5 As a consistency check we verify that these integrated mul- 
tipoles yield the same result as when calculating directly in the 
IHamiltonl J 1992T) algorithm. For discussion on binning issues with 
actual data see Appendix [El 



tegrating £(/i, s), and defining the radial £|| as /i > 0.5 and 
transverse £j_ as /i < 0.5. 

To calculate the "approximated" clustering wedges 
£\Vedge' that is, only taking into account the monopole 
and quadrupole contributions, we use Equation (|11|) . When 
adding the term of the hexadecapole fwedge' we ac ^ the 
contribution of Equation (|10[) . When using the formulae of 
iMatsubaral (|2004l ) (panel b) we also calculate Cwedge 6 ' 8 '' by 
following a similar procedure using ^6 and ^g. Results are 
given in terms of (^edge - Ccdg° XlmatCd )/°"«> where <r f cor- 
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Figure Bl. Joint 2a constraints on the H — D\ plane obtained 
when analyzing the AP effect for various clustering widths Afi 
in real-space (upper panel) and velocity-space (lower panel). For 
Afj, = 0.5 we test the effect of the correciton terms C\\ j_ (short- 
dashed blue; same results as in Figure [1J and without (thick dot- 
dashed black; Cm j_ = 0). (Ignoring intermixing terms means £| 
depends solely on H and £j_ on D\.) As indicated in the leg- 
end for each statistic, the symbols arc the most likely 2D values, 
and the crosses show the marginalized la results. We see that 
not including C|| t x biases the results, and underestimates uncer- 
tainties. We also show results for analyses of thinner wedges of 
Afi = 0.25, 0.125 (as indicated in legend). For these results we do 
not use correction terms in the analysis. We notice a clear trend 
of less bias as Afi decreases, which is expected as the intermixing 
terms are less important. The increase in uncertainty reflects the 
fact that we are using less information. In velocity-space this is 
seen very sharply by the elongation along the H direction. 



responds to the y/Cu obained from the 0.16 < z < 0.44 
LasDamas mocks. 

Due t o the fact that according to the prescription of 
Irlamiltonl (|l992l ) s) contains contributions only from 
even multipoles up to the hexadecapole, this £ = 4 con- 
tribution is critical for an absolute definition of the wedges. 
By neglecting this term the wedge approximations fledge 
are inaccurate at the ~ 1.5% level (in terms of fractions 
^approxy^exactj at barionic acoustic feature, which corre- 
sponds to 4% of <j£. 

For the formulae of iMatsub ara (200j) we follow a simi- 
lar procedure, where we also note contribution due to higher 



multipoles. We notice improvement from £ 



(0,2) 



to £ (0 ' 2 ' 4) 



and even further improvement with tiweAgf ^ ■ The ^6 



is 



contributions are not expected from the linear squashing ef- 
fect, but rather are a result of wide observer angle effects. 
We verify that ^6 and £s vanish when running the algorithm 
at z = 3. 

Comparin g the result s for obtained using the 

pr escription of Matsubara (|2004T ) to those from the recipe 
of lHamiltonl (|l992T l. we note that the former asymptots to 
zero fas ter. This can be e xplained by the fact that the £4 
term in iMatsubaral (|2004l ) is suppressed due to the rise of 
higher order terms. This is clearly seen in Figure [10] where 
£e is comparable to £4 at the baryonic acoustic feature and 
surpasses it at larger scales and £s is comparable to £4 at 
s > 130/i _1 Mpc. 

To conclude, these tests show that the clustering wedges 
can be accuratedly described in terms of the multipoles, and 
hence can be used as an alternative basis. 



D2 Analysis of mock catalogues 

Here we perform a similar analysis as in the previous section 
but using the results for the mean multipoles and wedges 
obtained from the 160 LasDamas mock catalogues described 
in £12.21 in real-space and velocity-space. As these include 
velocity-dispersion effects, as well as the expected wide-angle 
effects, we study even multipoles up to the eighth order. 
In practice, here all clustering multipoles and wedges are 
estimated by int e gratin g £(/x, s) which is calculated with the 
lLandv fc Szalavl (|l993h estimator. In Appendix[E]we discuss 
an alternative approach to calculating these ID statistics. 

We find that the monopole and quadrupole approxima- 
tions fledge appear to be sufficient to describe the clustering 
wedges on all scales in real-space. In velocity-space they ap- 
pear to be sufficient for scales s > 70/i _1 Mpc, as seen as 
thick short dashed lines in Figure IDlb .d. In real-space, of 
course, the only term is the monopole, but we do obtain 
higher order multipoles resulting from noise, and on large 
scales from angular effects. In velocity-space adding even 
multipoles I > 2 improves results at s < 70 h~ Mpc. These 
I > 2 terms appear due to velocity-disp ersion effects on lar ge 
scales and dominate at s > 30 /i -1 Mpc (|Scoccirnarrcll2004f l. 

In Figure fTOl we focus on the mock I = 4, 6, 8 multipoles, 
where uncertainties are for the mock mean of 160 realiza- 
tions . These a r e com pared to those analytical ly obtained 
from Hamilton] (|l992l ) (£4 only: dashed line) and Matsubara! 
(|2004F (solid lines). The mock £4 result appe ars to have 
struct ure similar to that expected according to IMatsubaral 
(|2004t l (that is, inclulding wide-angle effects). We might be 
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seeing a baryonic acoustic feature at ~ 100 /i~ 1 Mpc. Al- 
though the large uncertainties indicate low significance, we 
note that similar real-space analysis yields a negative £4 (but 
consistant with zero at a level of 2a), indicating that the 
velocity- space result is not due to angular effects only. Fig- 
ure 11 in lSamushia et aD (|201lf ) shows similar trends in the 
mock £4. 

We conclude that the two wide clustering wedges (A/x = 
0.5) are defined fairly well by the monopole and quadrupole 
in velocity-space (and monopole only in real-space), and 
hence may be used as a alternative basis to these multipoles 
to project most of the information contained in £(p, s). 



APPENDIX E: £ ESTIMATORS 

For our £ ana l ysis of the mocks we have used the 
lLandv fc Szalavl (|l993l ) estimator which has been shown 
in various studies as the lowest var iance known (e.g, 
iKerscher et al. I l2000l : iKazin et al.ll2010ah . We examine two 
methods in which the ID £t multipole and wedge statistics 
can be calculated: 

(i) "Integrated": calculating £(/z, s) in 2D bins and inte- 
grating using Equation [7] for the multipoles and Equation [6] 
for the wedges. 

(ii) "Direct" : calculating the statistics in ID bins directly 
in the pair-counter. 

For the "direct" method the ID estimators can be generi- 
cally written as: 



/^direct / \ 



DD 3M (s) + ARstat(s) - 2 • DR atat (s) 



RR 



projection type 



(s) 



(El) 



The DD, RR, DR are the usual normalized data-data, 
random-random, data-random pair counts. By normalized 
we refer to the fact that it is common to reduce shot-noise 
effects of the random points by assigning many more random 
points N n than data iV D : r = Nr/Nd > 10 This means, 
that when one counts number of data-data points Ndd the es- 
timator requires it to be normalized to the random-random 
iV rr by DD = Ndd ■ r 2 , while for data-random this is 
DR = iV dr • r. 

The "stat" subscript indicates the required weighting 
for each statistic. For the multipoles stat= £, meaning the 
weight is the £ Legendre polynomials Ve(l-i)- Technically this 
means when counting data-random pairs of the multipole 
A^ r (s), e.g, it is increased by: N dT + = Wd ■ w r ■ Ve(fJ-), if a 
pair is within range of [s — 0.5As, s + 0.5As] and has a [i 
within [/j,— 0.5A/1, fi+0.5A/j]. The weights w, could be due to 
incompleteness in uniformity (angular or radial), otherwise 
they are unity. For the multipoles the denominator term is 
that of the monopole, meaning -R-R mu iti P oios = RRo- 

The wedge estimators are similar in the nominator, 
but the denominator is different. For example, data-random 
pairs Nj^ dgc for clustering wedge /i m i n < /i < ^i m ax is in- 
creased by Wd ■ given a pair in the separation range 
[s — A/2s, s + As/2] and within the wedge fi range, and zero 
contribution otherwise. Because the nominator "weighting" 



6 For LasDamas catalogues we use r ~ 30, and for the Horizon 
Run mocks wc use r ~ 2. 



is that of a top hat, the denominator term i?-R wc dgc in the 
same as in the nominator. 

In this study we focused on two wide A/j, = 0.5 wedges, 
but this can, of course, be generalized to finer widths. One 
simple sanity check for A/i = 0.5 wedges is to verify that 
the average of both wedges yields the monopole. 

These two methods ("direct" and "integrated") should 
yield a similar result given RR is not a function of (i, but 
only of separation s. If this is the case, the "direct" method 
should be preferred because the "integrated" yields lower 
S /N due to binning e ffects . 

ISamushia et alj (|201ll ) point out, however, both wide 
angle effects as well as observer angle effects on the current 
SDSS-II geometry. When analyzing RR(fi, s) of the Las- 
Damas SDSS-II mocks, we find a ~ 10% of la difference 
between RR(1, s) to RR(Q, s) for the baryonic acoustic fea- 
ture scale, indicating the survey is wider than it is broader. 
In the BOSS 0.16 < z < 0.6 mocks we obtain similar con- 
clusions only this time a ~ 5% of la effect. In this case, by 
using the "direct" method (Equation IE1|I , one is introduc- 
ing an angular weighting of RR(fi, s) / RRo(s) such that the 
RR(/j,, s) is disregarded. In other words, Equation (|E1[1 is not 
equal to Equation (jHJ). For example, the "direct" multipoles 
actually take the form: 



&(-) 



DDtjs) + RRi(s) - 2DR e {s) 
RRo{s) 

2£ + l f +1 DD(fi,s) + RR(^,s) -2DR(fi,s) 



d/i- 



RR(fi,s) 



RR(/i,s) 
RR{s) ' 



(E2) 



In the left panel of Figure IE1I we show the difference 
between the two methods for clustering £0,2,4 and A/i = 0.5 
wedges. The difference between "direct" and "integrated" is 
given in units of for the SDSS-II mocks. 

We clearly see an over estimation of the "direct" results 
in most scales for the multipoles. At the baryonic acoustic 
feature scale this effect in units of la is ~ 5% for £0, ~ 7% for 
£2 and ~ 8% for £4, with the differences for the last two cases 
increasing at smaller scales. We note that the uncertainties 
of £2 and £4 are much larger than that of £0 making the 

ratios £f«ct ^integrated laj . ger than ^ q{ tfae monopole 

The A^i = 0.5 wedges, on the other hand, are shown to 
yield much smaller differences. This is probably due to the 
fact that the wide wedges do not correlate pairs at fi = 
with those at p, = 1. 

Another comparison is shown on the right panel of Fig- 
ure IE II which is similar to that in Figure [DU i. As described 
in Appendix [D] both test the relation between clustering 
wedges and multipoles (Equation I lip , where the right panel 
of Figure [Ell corresponds to the "direct" method, and Fig- 
ure [Dlji shows the "integrated" measurements. We clearly 
see that the "integrated" method performs better. 

Although the "integrated" estimator method appears to 
be closer to the natural definition of multipoles and wedges, 
there is a tradeoff of increase in uncertainties due to binning 
effects (one is integrating of more noisy DD). Also, by us- 
ing the "integrated" method explained her e, one might be 
addin g observer angle effects discussed in ISamushia et al.l 
|201l| y For this reason, we conclude that the user should be 
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Figure El. The left panel shows differences between the statistics obtained when using the "direct" method (£ dlroct ; ID binning in 
pair-counter) and the "integrated" method (£ mtc g ratod ■ integrating over 2D binning) in units of a^. Differences arrise due to wide angle 
effects for these clustering scales at (z) = 0.33. We clearly see that the wide wedges are the least affected. The right panel shows the 
relationship between the difference between the actual Afi = 0.5 clustering wedges to those statistics approximated when using multipolc 
combinations f M (Equation |6j in units of a^. In all cases shown here we use the "direct" measurements, while in Figure [DlH we show a 
simliar comparison for £ lnte g rated . \y e clearly see that the latter test is more successful, indicating the "integrated" method to be a more 
natural technique. 

aware of the two methods, and choose accordingly. We rec- 
ommend applying both and comparing the obtained results. 
The "direct" method might be preferable if the user ads the 
RR([i, s) I RRq(s) weighting into the multipole & model. For 
the A/i = 0.5 wedges we show no substantial difference in 
the actual £, so one would choose the one that yields lower 
uncertainties. 
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